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ABSTRACT 

The V-K colours along the minor axes of spiral galaxies typically change from red to 
blue with increasing distance, giving the impression that the near side is systematically 
screened by dust. Such a preferred orientation for dust screens is unlikely. Here we show 
that common extinction from the embedded dust layer in an exponential disk has the 
same effect, making the near side systematically redder as the inclination increases. 
The galaxy NGC 2841 is modelled as an example, where the V-K profile is profoundly 
asymmetric and actually step-like across the centre. We predict that the minor axis 
emission profile of the same dust in the far-infrared, at wavelength A ~ 200/im, will be 
much more symmetric than the optical profiles, implying nearly equal column densities 
of dust on both sides of the minor axis. 

Key words: galaxies: spiral - galaxies: ISM - galaxies: fundamental parameters- 
ISM: dust, extinction - infrared: galaxies - galaxies: individual: NGC 2841 
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1 INTRODUCTION 

The recent development of K'-band (2.1 pm ) and K- band 
(2.2 /xm) imagery for nearby galaxy disks led to the un- 
expected discovery of large V-K' and V-K colour gradients 
along the minor axes when the inclinations exceed ~ 60°. 
The "Evil Eye" Galaxy NGC 4826 exhibits a sudden increase 
by ~ 1 magnitude from V-K=4 mag. to V-K=5 mag. in the 
dust attenuated 'screen' (Block et al. 1994a). Similarly large 
V-K jumps are in other highly inclined spirals, such as NGC 
3521 (Panel 188 in Sandage & Bedke 1994). Explanations 
for this effect vary from an intervening dust screen (Witt 
et al. 1994) to high latitude scattering from dust (Block et 
al. 1996). The problem with the screen model is that the 
apparent dust is always on the near side of the galaxy, even 
though we expect random orientations for galaxies in space. 
The problem with the high latitude dust model is that there 
is no ubiquitous evidence for such dust, although extrapla- 
nar dust in specific examples has been detected (eg. Howk 
& Savage 1997). 

This paper shows how absorption and scattering at V 
and K bands by dust embedded in an exponential disk can 
produce observed colour gradients which may be remark- 
ably asymmetric or even step-like. The reason is that the 
near side of an exponential disk is brighter at smaller ra- 
dius behind the dust in the midplane, while the far side is 
brighter at smaller radius in front of the dust in the mid- 
plane. Thus the near side has a higher fraction of the same 
total light blocked by dust. 



Other models have considered similar effects. Kodaira 
& Ohta (1994) and Ohta & Kodaira (1995) measured 
and modelled V and J band differential extinction profiles 
(=— ln[Ifar/Inear]) along the minor axes in several galax- 
ies, without including scattered light; they did not discuss 
colour gradients specifically, although the existence of such 
gradients can be inferred from their analyses. Byun, Free- 
man & Kylafis (1994) considered B and I band minor axis 
profiles with scattering included, but also did not specifi- 
cally address the resulting near-far colour differentials (al- 
though again they can be inferred from their separate B 
and I profiles). Kuchinski & Terndrup (1996) plotted J-K 
colour profiles along the minor axes of several galaxies, but 
discussed and modelled only the spherical bulge systems, in- 
cluding scattered light. Here we determine the line-of-sight 
V-K colour gradient in disk+bulge galaxies using radiative 
transfer with direct and single-scattered light. 



2 MODELS 

Stellar and gaseous disks are modelled with exponential 
brightness and density profiles in both the radial and per- 
pendicular directions, following the work of others referenced 
above. For the stellar volume emissivity, we use 



j(r, z) — exp —r 



\ z , 

~) + (i-RyR 2 b ) 15 



G) 



and for the dust extinction (scattering + absorption), 
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K\(r,z) = K X:D (0,0)exp -r 



M,b(0,0) 



(2) 



All distances are normalized to the disk scale length; r is 
the distance to a point measured parallel to the disk, R is 
the total galactocentric distance, in three dimensions. The 
volume emissivity is normalized to unity for each wavelength 
because we are concerned primarily with differential (near- 
far side) colour gradients. The intrinsic stellar colours do not 
matter as long as they are the same at each radius on the 
near and far side of the disk. Radial colour gradients in the 
disk do not matter for a near-far comparison either. 

The equation of radiative transfer gives the intensity I 
measured by an observer outside the galaxy at wavelength 
A: 



I\{i)= I dx(j [r, z] + j scat [r, z] ) 

JXQ 



x exp 



K\[r, z]dx' 



(3) 



The distance along the line of sight is x, increasing away 
from the observer; xq and x\ are the near and far limits of the 
integral through the galaxy, taken here to be the hypotenuse 
of a triangle with a maximum total in-plane radial extent of 
5 scale lengths and a maximum total half-thickness of 3 scale 
lengths. The argument of I\, i, is the galaxy inclination. 
The volume emissivity from single scattering is 



jscat(r,z) = K X (r, z)A x 



dl 



r/2 



dfecos b(f>{6) 



r/2 



/ j(r, z)d£exp - / K\[r, z]d£' 
V Jo 



(4) 



where X2 is the maximum distance along a path through 
the galaxy, I and b are the galactic longitude and latitude 
as seen by an observer at point x in the integral for I\. is 
the scattering phase function (Henyey & Greenstein 1941), 

1 1 „2 



4n (1 + 9 2 - 2.gcos<9) 1 - 5 
and is the scattering angle, given by 
cos 6 — — sin b cos i + cos b sin i cos I. 



(6) 



In equation (4), A\ is the albedo at wavelength A. While 
most interstellar grains were thought to be substantially 
smaller than near-infrared wavelengths (Mathis, Rumple & 
Nordsieck 1977), the upper size limit of dust grains - to 
which the near-infrared albedo is extremely sensitive - has 
been revised (see Kim, Martin & Hendry 1994). The first ex- 
tragalactic determination of the near-infrared dust albedo 
was made by Witt et al. (1994), who concluded that the 
albedo of dust grains at V and K may be identical. This high 
near-infrared albedo is alluded to by a number of diverse 
observations (Lehtinen & Mattila 1996 for the Thumbprint 
Nebula; Block 1996 for the Whirlpool Galaxy M51 and its 
companion; Pendleton, Tielens & Werner 1990 and Sellgren, 
Werner & Dinerstein 1992 for infrared and classical reflec- 
tion nebulae). We choose values for the optical and near- 
infrared albedos of dust grains to both be A — 0.6 in the V 
and K passbands. 

Furthermore, dust grains are predominantly forward 
scattering in the optical (eg. Whittet 1992) while near- 
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Figure 1. (top) V-K colour excess profiles along the minor axes of 
model galaxies with inclinations of 50°, 60°, 70°, and 80° (E[V - 
K] increases with inclination). Solid lines include scattered light. 
This case has no dust distributed in the bulge, (bottom) V-K 
colour excess differences between the near and far sides shown as 
functions of minor axis distance. 



infrared values of g at K lie between 0.1 and 0.6 (Witt et 
al. 1996). We choose g — 0.5 in V and K. There is only 
a weak dependence on g as far as albedo determinations 
are concerned (Fig 2 in Witt et al. 1996), and our models 
with different values of A did not change much; lower A 
gave about the same results as slightly higher extinctions, 
approaching the results with no scattering, as shown in the 
figures by dashed lines. 

The intensity 7a was determined from these equations 
as a function of distance from the galactic center for various 
galaxy inclinations. For each distance and inclination, the 
reddening difference between the near and far sides of the 
disk was evaluated in magnitudes: 

„ - E(V - JOfa* = (7) 

2.51og 10 (^) . (8) 

\lK/ far 

For most of the models shown here, the perpendic- 
ular scale heights are taken to be z s = 0.1 and z s — 
0.2, and the bulge scale length is Rt = 0.1 (Kodaira, & 
Ohta 1994). The relative volume emissivity of the bulge is 
e b ~ 0.5 (Lb/Ld) Zs/Rl (Kodaira & Ohta 1994). The ra- 



E(V - K) t 
-2.51og 10 (^) 
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Figure 2. Same models as in the previous figure, but with dust 
distributed throughout the bulge region, in three-dimensions. 



tio of bulge to disk luminosity is Lb/Ld = 0.1 at K band 
(Byun, Freeman & Kylafis 1994) and 0.1 x i(r - 4AM at 
V band for colour difference A M = 0.5 mag between the 
bulge and the disk. The central extinctions for the disk 
are taken to be kv,.d(0, 0) = 10 or 15 optical depths (o.d.) 
kpc -1 for the models shown by the figures, and for the bulge 
are kv,b(0,0) = 0, or 10 o.d. kpc -1 (1 optical depth = 
0.4/ log e = 0.92 magnitude.) In all cases, kk = O.lftv. The 
value of k,v,d{0,0) = 10 o.d. kpc -1 might be appropriate 
for the Milky Way: it gives a V-band extinction of 1.25 mag 
kpc -1 at r — 2 disk scale lengths, which is about the ex- 
tinction at the solar radius. Many other values for these 
quantities were modelled too, and some of the results will 
be discussed for comparison. 

Figure |l|(top) shows the intensity versus distance along 
the minor axis for a model with a central disk extinction 
of 10 o.d. kpc -1 in V and no additional extinction in the 
bulge. The projected distance along the minor axis, xcosi, 
is on the abscissa, in units of the intrinsic disk scale length; 
negative values of this distance represent the near side of the 
galaxy. Both the bulge and disk emissions are included here, 
but there is no bulge extinction in this model. Cases with 
scattered light are shown by solid lines, and cases without 
scattered light are shown by dashed lines. Inclinations of 
50°, 60°, 70°, amd 80° are plotted; the near-side reddening 
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Figure 3. Minor axis V-K' profile for the galaxy NGC 2841, 
showing redder colours on the near side (left) than the far side. 
The axis for the colour difference is on the left. Also shown are 
the V and K' intensity scans across the minor axis, using the axis 
on the right. The dashed line is K'. 



increases with inclination. Any intrinsic disk colour gradient 
that might be present in the stars is ignored. 

The V-K colour differences between the near and the 
far sides are shown in the bottom panel of figure There is 
a maximum in differential reddening at about 0.2 disk scale 
length, where the visual extinction is ~ 8 o.d. kpc - . Intrin- 
sic disk colour gradients would not matter for this diagram, 
as long as the near and far sides of the disk have the same 
intrinsic colours at the same radii. 

Figure [] shows the same model as figure |l| but with ex- 
tinction in the bulge. For the assumed central bulge extinc- 
tion of k,v,b(0, 0) = 10 o.d. kpc -1 , the extinction at a radius 
of 1 disk scale length is 10 o.d. kpc -1 (l + [1/0. if) ~ 3/2 = 
0.01 mag kpc -1 , which is much smaller than the disk extinc- 
tion at this radius. Nevertheless, dust in the bulge clearly 
reddens the far side of the V — K colour profile, making the 
V-K colour difference at the bottom of the figure about a 
factor of 2 less than in figure ^. 

Figure ^ shows the observed V-K' scan along the minor 
axis of the inclined flocculent galaxy NGC 2841 (studied in 
more detail by Block et al. 1996). The strong colour gradient 
between the near (left in the figure) and far sides is evident, 
as discussed in our previous paper. The models in figures 
[l] and ^ do not reproduce this gradient well because NGC 
2841 has a central hole in molecular gas inside 1 disk scale 
length (Young & Scoville 1982), and it has a low HI column 
density there too. Thus there is probably a dust hole inside 
1 scale length. This hole was also evident in a scan of V-K 
colour versus distance along the major axis of NGC 2841 
(Block et al. 1996), which shows a dip toward bluer colour 
inside 1 disk scale length. This blueness approaches the true 
colour of the inner disk, without extinction. In addition to 
the dust hole, NGC 2841 also has a relatively high extinction 
perpendicular to the disk compared to the Milky Way (Block 
et al. 1996), so we need to consider slightly larger k than in 
the previous figures. 
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Figure 4. Model appropriate for NGC 2841 with four inclinations 
i as before. In this case the inclination i = 70° , which is the second 
from the top, most closely matches the inclination of NGC 2841. 
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Figure 5. Model with a disk dust distribution as in Fig. 4, but 
with dust also in the bulge, off the plane. This case is not a good 
match for NGC 2841 because the bulge dust reddens the far side 
of the disk, and produces too high a colour differential there. 



Figure |4| shows the minor axis profile from a model with 
1.5 times the extinction as in figure |l| but with a central 
hole made by setting the dust density profile equal to re~ r 
instead of e~ r for radius r, measured in units of the scale 
length. This function gives a peak in the gas density at r = 1 
with a hole inside this and a nearly exponential disk beyond. 
The model has the same four inclination angles as in the pre- 
vious plots, for general use in other studies; the inclination 
of NGC 2841 is ~ 68°, which is essentially the same as for 
the third curve up from the bottom (which is for 70°). With 
a hole in the model and no dust in the bulge, the fit to the 
data is good, i.e., the near side is redder than the far side by 
~ 0.4 mag., and the far side colour excess distribution is rel- 
atively fiat. The small scale variations in the data for NGC 
2841 are presumably from dust lanes, which are not part of 
the present model (but see Block et al. 1996 for radiative 
transfer fits to these dust lanes). 

Figure [| shows the same disk-hole model as in figure 
but with extinction in the bulge as in figure ^| This model 
does not fit the data for NGC 2841 because the far side 
is too red with the additional extinction from the bulge, 
giving the far-side part of the curve on the top of figure |fj| a 
downward slope, unlike the observations of NGC 2841 and 
the corresponding curve in figure ^, where this far-side E(V- 



K) scan is flat. The dusty-bulge model also has an increased 
reddening in the central region, which becomes even more 
prominent for face-on orientations. Any larger Kv,b(0,0) 
makes this central reddening unacceptably large even for the 
inclination of NGC 2841. Thus the upper limit to the amount 
of bulge dust in NGC 2841 is about the value in figure 
which corresponds to less than 0.01 o.d. kpc - at one scale 
length, or an equivalent hydrogen density of < 0.007 cm -3 , 
using the standard conversion of dust to gas (Bohlin, Savage 
& Drake 1978). Evidently our previous suggestion (Block et 
al. 1996) that the near-far colour difference in NGC 2841 is 
the result of bulge scattering by dust high off the plane is 
not correct. 

The variable geometry of dust and starlight along the 
line of sight makes the relationship between color excess and 
extinction non-linear, unlike the common interstellar case 
where an absorbing cloud is in front of a single star. Fig- 
ure ^ plots the reddening versus the optical depth through 
the disk for the Milky Way and NGC 2841 models with no 
bulge dust (from Figs. |l| and ^) and with scattered light 
included. The projected radius along the minor axis varies 
along each curve, from the near side of the disk (solid line) 
to the far side (dashed line). The same four inclination an- 
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Figure 6. Color excess versus optical depth through the disk for 
the Milky Way and NGC 2841 models shown previously in figures 
1 and 4. The lack of a simple linear relationship results from the 
varying distribution of dust and starlight on the line of sight. 

gles are shown. The relationship between reddening and ex- 
tinction is double-valued because the near side has more 
reddening per unit dust column density than the far side. 



3 DISCUSSION 

Differential reddening between the near and far sides of a 
galaxy disk may be understood as follows: The radial light 
gradients in the bulge and disk cause the near side to have a 
brighter light source behind most of the disk dust, while the 
same gradients on the far side cause it to have its brighter 
light source in front of the disk dust. Thus the disk dust has 
more light to block on the near side than the far side. Con- 
sequently the near side is dimmer (Kodaira & Ohta 1994; 
Byun, Freeman, & Kylafis 1994) and significantly redder. 
Models with the same scale height for the dust and stars 
(not shown) have even stronger colour differentials because 
more of the starlight is extinguished. Models with more disk 
extinction also have more differential reddening, and mod- 
els with dust in the bulge have less because the bulge dust 
makes the far side redder. 

Minor axis colour gradients are present in many of the 
galaxies in Wray's Color Atlas (1988) and are even more 
prominent in V — K images of inclined galaxies (e.g., see 
Peletier and Balcells 1997; Thornley 1996 & 1997; Grosb0l 



— i 1 r 

NGC 134 (from Alton et al. 1998) 




-100 -50 50 100 

Radius (arcsec) 

Figure 7. Surface brightness along the minor axis of the galaxy 
NGC 134 at 200 urn, from Alton et al. (1998), showing symmetry 
in cold dust emission even though this same dust presents a strik- 
ing asymmetry with foreground "screening" in optical images. 

& Patsis 1998; Block et al. 1994a). Thornley (1996) found 
that in NGC 5055, optical depths suddenly jump from 0.5 
on the far side, to 4-5 on the near side. 

The models presented here show that such colour gra- 
dients are a natural result of extinction from well mixed 
dust embedded in the disk, and it is predicted that for large 
enough inclinations, particularly in galaxies with inner-disk 
dust holes, the near-far asymmetry in V-K profiles may even 
be step like. Near-side dust screens are not needed, even for 
the dramatic jump of E(V - K) ~ 1 mag in NGC 4826. 
There is some sensitivity of the colour gradient to the pre- 
cise distribution of dust, such as the presence or lack of a 
central hole, and this may be useful in modelling disks. 

Another point to stress is that the far side, although 
bluer, cannot be considered dust free or unobscured - terms 
that continue to be used in the literature. For example, the 
Andromeda spiral M31 shows prominent dust lanes on the 
far side as well as the near side in both the IRAS 60 /im and 
100/xm maps (eg. Habing et al. 1984; Walterbos & Schwering 
1987). NGC 2841 also has arms of dust in K' images on both 
sides (Fig lc in Block et al. 1996 and plate 3 in Block 1996), 
as does M51 at 15/^m (Sauvage et al. 1996, Block et al. 1997). 
It is simply that when the disk is well inclined to the line 
of sight, the near-side dust appears to be much thicker than 
the far side dust, even when both the dust and the stellar 
column densities are the same on each side. 

These results lead to a prediction about the emission 
from cold (T~ 20K) galactic dust: the far-infrared minor 
axis luminosity profile across a galaxy, at wavelengths of 
~ 180 /im or longer, should be symmetric, even if the optical 
profile is so asymmetic that it appears to be screened by dust 
on one side. That is, optically screened galaxies will be found 
to have nearly equal column densities of cold dust on both 
sides of the minor axis. 

A good example is NGC 134, which was cited by 
Sandage & Bedke (1994) as displaying a "particularly 
strong" dust asymmetry between the near and far sides (see 
Panel 194 in their Atlas). This galaxy has just been ob- 
served at 200/im (Alton et al. 1998) by the Infrared Space 
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Observatory. The far-infrared emission profile of NGC 134, 
shown in figure Q is in fact almost perfectly symmetric 
across the minor axis (Alton, private communication) . Thus 
equal amounts of cold dust, which accounts for ~ 90% of 
the dust mass (Block 1996; Block et al. 1994a), reside on 
either side of the minor axis, despite its apparently striking 
foreground dust screen. 

Another implication of the model results, illustrated 
by figure ^| is that galactic opacity is not simply related 
to color excess. The conversion from color excess to opac- 
ity depends on inclination and position in the disk. For 
galaxies well inclined to the line of sight, there are many 
E(V-K) values for the same optical depth, with extrema 
on the minor axes and intermediate values elsewhere. If 
this is not taken into account, one might erroneously find, 
on the basis of E(V-K) colour excesses, that the bluer far 
side contains less dust than the redder near side - and 
will incorrectly compute differences in the dust masses of 
the near and far sides - although both sides contain equal 
amounts. Figure ^ does suggest, however, that the average 
of the color excesses for the near and far sides might have 
some simple relation to the opacity at that radius when 
there is no inner hole. From the bottom panel of this fig- 
ure, we obtain r ~ 2.5 {E[V - K] near + E[V - K] far ) /2. 
Converting this to visual magnitudes gives Av ~ 
2.7 (E[V - K] near + E[V - K] far ) /2. Note that this con- 
version factor is much larger than the conventional factor 
from V-K color excess to Ay, based on foreground extinc- 
tion for single stars, which gives Av ~ 1.09E(V —K) (Savage 

6 Mathis 1979; Whittet 1988). The larger factor here arises 
because nearly half of the dust on the line of sight is on 
the far side of the disk, and this dust occults relatively few 
galactic stars. 

Our results might also have a practical application. We 
have noted that when a minor axis colour gradient is found, 
the red side should be nearer. Thus the sense of dominant 
rotation in a disk galaxy can be determined from the colour 
gradient if trailing spiral arms are present. Such determina- 
tions may be useful for crowded fields of distant galaxies in 
an attempt to look for alignment or other attributes of spin 
orientations without the tedious job of taking slit spectra to 
get rotation curves. This characteristic of disk colour might 
also be useful in searching for leading spiral arms in galaxies 
with rotation curves. Such arms apparently occur in some 
accreting galaxies (NGC 4622: Buta et al 1992, see also Fig 

7 in Block et al. 1994b and plate 4 in Bertin & Lin 1996; 
NGC 4826: van Driel & Buta 1993). 
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